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New measurement by CELSIUS- WASA Collaboration on the pp — > pnir + reaction re- 
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veals clear evidence for the presence of the Roper resonance A* (1440) which has been 
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ignored in previous theoretical calculations. In this article, based on an effective La- 
grangian approach and available knowledge on the Roper resonance, we investigate the 
role of the Roper resonance for the pp — ► pnir + reaction. It is found that the contribu- 
tion from the Roper resonance A* (1440) becomes significant for kinetic energy above 1.1 
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GeV, consistent with the new experimental observation. The t-channel a-meson exchange 
is dominant for the production of the Roper resonance. 

PACS: 13.75. Cs; 14.20. Gk; 13.30.Eg 

keywords: pp collision; Pion production; iV*(1440) resonance; Effective Lagrangian ap- 
proach 

1 Introduction 

Ever since the first even-parity excited state of the nucleon, the Roper resonance iV*(1440) 
was first deduced by nN phase shift analysis, its structure has been arousing people's interests 
intensely all the time — it is lighter than the first odd-parity nucleon excitation, the S , n(1535), 
and also has a significant branching ratio into two pions. Up to now, although the existence of 
the Roper resonance is well established (4-star ranking in the particle data book), its properties, 
such as mass, width and decay branching ratios etc., still suffer large experiment uncertain- 
ties [1]. In classical quark models, the Roper resonance has been associated with the first 
spin-parity J p = l/2 + radial excited state of the nucleon [2, 3, 4, 5]. In the bag model [6] and 
in the Skyrme model [7], a surface oscillation, also called breathing mode, has been predicted 
to interpret the Roper resonance as a monopole excitation of the nucleon. Furthermore, it has 
also been supposed to relate to a hybrid nature which means a gluonic excitation state of the 
nucleon [8, 9, 10]. The recent theoretical works [11, 12] found that the Roper resonance was 
dynamically generated from the meson-nucleon interactions. Nevertheless, the model predic- 
tions always reach either larger value for the Roper resonance mass or much smaller one for its 
width and also meet difficulties to explain its electromagnetic coupling [13]. 
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In the early years, our knowledge on Roper resonance was mainly based on ttN and 
7iV experiments, however, it is excited so weakly thus buried by the strong A peak in these 
experiments. Another difficulty in extracting the Roper resonance information from these 
experiments is the isospin decomposition of 1/2 and 3/2 [14]. In Refs. [15, 16], it was pointed 
out that the decays of J/ip — > NNtc and J ftp — > NNtctc provide an ideal place for studying 
the properties of nucleon resonances, since in these processes the nN and 7T7riV systems are 
limited to be pure isospin- \ due to isospin conservation. The result from BES Collaboration on 
J/ip — > pn~n+ c.c. reaction showed a clear N* peak in the Nn invariant mass spectrum around 
1360 MeV/c 2 which gave the first direct observation of the Roper resonance peak in the nN 
invariant mass spectrum [17]. Nevertheless, since the Roper resonance has a strong coupling 
to aN, studying the Roper resonance from other production processes, i.e., aN — > iV*(1440), 
is necessary, as also suggested in Refs. [18, 19]. 

Recently, the CELSIUS- WASA Collaboration performed a measurement on the pp — > 
pnn + reaction with proton beam of kinetic energy T P = 1.3 GeV. They also observed a clear 
resonance structure in the invariant nir + mass spectrum for the Roper resonance atMsi 1360 
MeV/c 2 with a width of about 150 MeV [20, 21]. These values agree well with the BES result [17] 
and the pole position of the Roper resonance from nN phase shift analyses [1]. The pp — > pnn + 
reaction has the largest cross section for pp collision in the intermediate energy region around 
T P = 1.3 GeV. It can be further studied by WASA-at-COSY and HIRFL-CSR facilities with 
much higher precision and statistics. This reaction opens a new window for studying the Roper 
resonance. 

However, up to now, theoretical study on this channel in the intermediate energy region 
is scarce. In Ref. [22], Engel et al. performed a fully covariant calculation for the total cross 
section of the pp — > pn-K + reaction. However, unaware of the significant role of the Roper 
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resonance for this reaction, their calculations only included the A resonance and the off-shell 
nucleon pole. Furthermore, in their calculations they did not consider the final-state interaction 
of pn. In Ref. [23], the authors performed their calculations for T P = 0.8 GeV by including 
only the A resonance. 

Inspired by the new observation of the Roper resonance by the CELSIUS- WASA Collabo- 
ration [20, 21], here we investigate the role of the Roper resonance for the pp — > pnn + reaction, 
based on an effective Lagrangian approach and available knowledge on the Roper resonance. 

In next section, we will give the formalism and ingredients for our calculation. Then 
numerical results and discussion are given in Sect. 3. 

2 Formalism and ingredients 

We study the pp — > pnn + reaction in an effective Lagrangian approach. All the basic Feyn- 
man diagrams involved for the reaction are depicted in Fig. 1. The so-called "pre-emission" 
diagrams [22] are not shown. But their contributions are included in actual calculations. In ad- 
dition to the A(1232) and the off-shell nucleon pole included in previous works [22, 23, 24, 25], 
the Roper resonance is added in our calculation. For the production of the Roper resonance and 
the intermediate virtual nucleon, the n° and a-meson exchanges are included while the p-meson 
exchange is found to give negligible contribution and is dropped. For the production of the 
A (1232) resonance, we only include 7r-meson exchange since iso-scalar meson exchanges cannot 
contribute and according to Ref. [23] any attempt to include the p-meson exchange worsens the 
agreement with experiments of pp — > nA ++ reactions. A Lorentz covariant orbital-spin(L-S) 
scheme [26] are used for the effective AirN, N*aN, and N*tcN vertices. 
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2.1 Meson-Baryon-Baryon (Resonances) vertices 

The effective Lagrangian densities involved for describing the irNN and aNN vertices are, 

C-nNN = -^^-un^^t ■ d^^iiN, (1) 

£<tNN = daNNUN^PaUN- (2) 

At each vertex a relevant off-shell form factor is used. In our computation, we take the 
same form factors that used in the well-known Bonn potential model [27]: 

A 2 2 

TpNN( h 2 \ _ A M ~ m M (n\ 
r M K^M) — A 2 _ 1,2 ' \°) 
A M K M 

where &m, n^M and Am are the 4- momenta, mass and cut-off parameter for the exchanged 
meson (M), respectively. In the Bonn meson-exchange model [27] for the nucleon-nucleon 
interaction, k° M = 0, hence k 2 M = — \k M \ 2 is negative. For the pp — > n + pn, it can easily verify 
that for the t-channel meson exchange, the k 2 M still keeps negative. Therefore the off-shell 
form factors here always give a suppression factor. The coupling constants and the cut-off 
parameters are taken as the following ones [27, 28, 29, 30, 31, 32]: g 2 NN /^ = 14.4, A„- = 1.3 
GeV, g 2 . NN /Arr = 5.69 and A CT = 1.7 GeV. Note that the constant g n NN is related to f w NN of 
Eq. (1) as QttNN = {UNN/m^)2m N [33]. 

To calculate the amplitudes of diagrams in Fig. 1 with the resonance model, we also 
need to know interaction vertices involving A(1232) and A r *(1440) resonances. In Ref. [26], a 
Lorentz covariant orbital-spin scheme for N*NM couplings has been described in detail, which 
can be easily extended to describe the A(1232)Ntt, N*(1440)A% and A^*(1440)A^a couplings 
that appear in the Feynman diagrams depicted in Fig. 1. By using that scheme, we can easily 
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obtain the effective A(1232)A%,iV*(1440)iV7r, and N*(lU0)Na couplings: 

^7rWA(1232) = 9A(1232)NTrUNd tJ, tp 7r UA(1232)^ + h.C., (4) 
AW* (1440) = 9N* (1440)^^757^^^^* (1440) + h.C. , (5) 
ArAW*(1440) — 9N* (U40)NaUNlpaUN* (1440) + h.C, (6) 

with un, Wjv*(i44o), ma>, ipw, and ip a as the wave functions for the nucleon, A/"* (1440) resonance, 
A* (1232) resonance, n and cr-meson, respectively. 

For the relevant vertices involving A(1232) and A/"* (1440) resonances, the off-shell form 
factors are adopted as below: 

FM{k 2 M ) = ( %Z™h n > (7) 
where n=l for the Roper resonance and n=2 for the A(1232) resonance [18, 22, 34], with 
Am being the corresponding cut-off parameters which are: = 0.59 GeV, A^ ( - 1440 ' ) = 1.3 
GeV [18] and ( - 1440 - ) = 1.5 GeV [34]. Here a quite severe cutoff A^ is needed to reproduce 
data, in consistent with the previous study [23] of pp — > nA ++ reactions. 

It is worth noting that we also take the virtuality of the intermediate nucleon into account 
by introducing a form factor which is taken as in Refs. [35, 36, 37] 

with Aiv= 0.6 GeV. 

2.2 Coupling constants for the intermediate resonances 

The A(1232)A% and A r *(1440)A% coupling constants are determined from the experimentally 
observed partial decay widths of the A and Roper resonances. The general formula for the 
partial decay width of A(1232) or V*(1440) resonance decaying into a nucleon and a pion is as 



the following 



d 3 p N m,N d 3 p n 1 



iv = lM ^ M V (qR - PN - P ^^^^-, (9) 



where Mr^ntt stands for the total amplitude of the initial resonance decaying to a nucleon 
and a pion, the qn, p^ and p w are the 4- momenta of the three particles, E N and E n are the 
corresponding energies. With the effective Lagrangians described by Eq. (4) and Eq. (5), the 
partial decay widths can be calculated 

3g 2 N * (imN7T P C N ml(E N -m N ) 2 

1 TV* (i440)-JVir = — [ —j 1" 1 \Vn ) \ ( iU J 



1- M, 



AT* (1440) 



with 



p5T = 



\ 



( M /V*(1440) - ( m N + m ^) 2 ){ M N*(lAm) - ( m N ~ m ^) 2 ) 



AM 2 

N* (1440) 



(11) 



En = ^{pT) 2 + m 2 N , (12) 

r _ 9 2 A *(l232)NTr( m N + En)^) 3 

1 A*(i232)Af7r ToTnT ' { L6 > 

127rM A *(i232) 

with p™ here the momentum of the nucleon in the rest frame of A(1232). 

For the iV*(1440)iV<7 coupling constant, we get it from the partial decay width of iV*(1440) - 
Na — > N-n-n, which is given by 

, r YU 12 1711 1 1 (O \4 r4/ \ 

(14) 

with 

Mn'^No^N™ = ^g^a9N*NaF^* N (kl)u N (p U Si)— 5-— — p 2 • P3^AT* (p/V* , S N .) , (15) 

ft* - + im a Y a 

where M.n*->n<7->Nwx represents the total amplitude of iV*(1440) — > Na — > A^7r7r. p/v* and 
denote the 4-momenta of the iV*(1440) resonance and the intermediate a meson; m a and Y a 



are the mass and full width of the a meson, which we take m a = 550 MeV and r CT = 500 
MeV; pi, m 1 , and E 1 stand for the 4- momenta, mass, and energy of the nucleon; Si and s N * 
the spin projection of the nucleon and the iV*(1440) resonance; p 2 , P3, and E 2 , E 3 stand for 
the 4-momenta and energy of the final two pions, respectively To calculate the amplitude of 
cr — > 7r7r appearing in the total amplitude calculation, we use the following type of vertex 

= fiW^VV • d"^M\r, (16) 

where the coupling constant g n7ra can be determined from the partial decay width T a ^ n7T 

jt- 1 TTa P?d 3 p 1 1 d 3 p 2 1 N4r4 , s „„s 

rfr ™ = ^'^— 1 (^^(^2^ (27r) 5 b*-n-P*\ (17) 

where M. g ^-k-k represents the total amplitude of o — > nn. pi, p 2 , and E±, E 2 stand for the 
4-momenta and energy of the final two pions, respectively 

With the experimental branching ratios [1] and the above formulae, we obtain all the 
corresponding coupling constants as summarized in Table 1. 

2.3 Propagators 

In our amplitude calculation of Feynman diagrams in Fig. 1, we also need propagators for the 
7r and cr-meson, and the intermediate A(1232) and iV*(1440) resonances with half-integer spin 
as well. For the n and a-meson, the propagators are 

K ir/a m 7T/a 

where k n / a are the 4-momenta of n and cr-meson, respectively; m 7T / CT are their corresponding 
masses. 

The propagators for the A(1232) and iV*(1440) resonances can be constructed with their 
projection operators and the corresponding Breit-Wigner factor [38]. The spin-(n+|) resonance 
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propagator can be written as 

[q) <p-M*+iM R r R ' li9j 

where l/(q 2 — M R + iM R T R ) is the standard Breit-Wigner factor; M R) q and T R are the mass, 
four momentum and full width of the resonance, respectively. p(™ + 3) is the projection operator 



p,i (9 ) = (#+Mfi) 



2M 



R 



(21) 



With the projection operators, the propagators for the intermediate A(1232)(| + ) and 
iV*(1440)(| + ) resonances can be easily obtained as the following 

qN*(U40)^ = j + M jV*(1440) ^ 

q 2 — Mjj*^ im ^ + lM N * ( 1440 ) Y N * (1440) 



/yA(1232)/ \ _ Py-Aq) 

^ w g 2 -Mi {1232) +^M A(12 3 2) r A(1232 ) 



(23) 



with 



PM = ( + M A(1232 ) ) \g - \ ltllv - lv% - — Y q^q v \. (24) 

3 3M A(1232 ) 3M A(1232) 

For simplicity we use constant width in the Breit-Wigner (BW) formula. The mass and 
width appearing in such constant-width BW formula are very close to their corresponding pole 
positions. In some detailed data fitting works, various energy dependent widthes are used in the 
BW formula. The mass and width appearing in such BW formula could then be very different 
from their corresponding pole positions. Usually the pole positions are better determined 
than the model dependent BW mass and width. For example, the BW mass and width given 
by PDG [1] with energy dependent width are (1420 ~ 1470) MeV and (200 ~ 450) MeV, 
respectively; while the corresponding pole positions listed in PDG [1] are (1350 ~ 1380) MeV 



and (160 ~ 220) MeV, respectively. Hence, for parameters appearing in the constant-width 
BW propagators, we adopt M A (i232) = 1-21 GeV, r A (i232) = 100 MeV, Mjv*(i44o) = 1.365 GeV 
and rjv.( 1440 ) = 190 MeV, as given by PDG for their corresponding pole positions [1]. 
As for the intermediate virtual nucleon, we use 



as its propagator, with M N and q being the mass and 4-momenta of the virtual nucleon, 
respectively. 

2.4 Amplitude and total cross section for pp — > pnn + reaction 

In order to show the structure of total amplitude in terms of the meson exchange and the 
intermediate resonances, we write down explicitly the total amplitude in our calculation for the 
pp — > pnn + reaction as the sum of sub-amplitudes 



M = A4(A(1232),7T + )+A4(A(1232),7r )+7W(Ar*(1440),7r )+7W(iV*(1440), ( 7) + 



where on the right-hand side of the equation, the resonances and mesons exchanged in the in- 
termediate states are written inside the bracket explicitly. Each sub-amplitude can be obtained 
straightforwardly with effective couplings and propagators given in former sections by following 
the Feynman rules. The explicit expressions of the amplitudes are given in Appendix. 

For the fmal-state-interaction(FSI) between the three outgoing particles, we only need to 
consider p-n FSI enhancement factor near threshold, since the interaction between nn + and pn + 
are dominated by the s-channel intermediate resonances. To study possible influence from the 
p-n final state interaction, we include it in our calculation by factorizing the reaction amplitude 



G N (q) 



4+ m N 

q 2 - m 2 N 



(25) 



M(p(938), 7T°) + M(p(938), a), 



(26) 
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as 



A = M(pp^pnir + )T pn , (27) 

where M.{pp — > pmr + ) is the primary production amplitude as discussed above, T pn describes 
the p-n final state interaction, which goes to unity in the limit of no FSI. The enhancement 
factor T pn is taken into account by means of the general framework based on the Jost function 
formalism [39, 40] with 

T m = £±£ (28) 
k — la 

where k is the internal momentum of p-n subsystem, and the a and (3 are associated with the 
effective-range parameters via [41] 



1 / 2r 1 / 2r 

a=-l-Jl--, /? = -(! + A/1 )■ (29) 

r V a r V a 

According to the FSI analysis of Ref. [41] on the reaction pp — > pnir + at the beam energy 
of T p = 0.95 GeV, the triplet 3 Si np state dominates the np invariant mass spectrum at the 
low mass end with negligible singlet 1 S'o state contribution. So in the present work, we assume 
the total amplitude with triplet 3 Si np state dominating the np near threshold region, and 
multiply the total amplitude with the FSI factor of triplet 3 S'i np state with the parameters 
adopted from Ref. [41], i.e., 

at = 5.424 fm, r t = 1.759 fm. (30) 

Then the calculation of the invariant amplitude square \A\ 2 and the cross section a{pp — > 
pnn + ) are straightforward, 

, / +x l m pv^v^ A i2 m pd 3 P3 d 3 p w m n d 3 p n 4 

da{pp^pn^) = -— 2^2^\A\ — — — 5 { Pl + p 2 - p 3 - p n - p n ), (31) 

s t s f 3 w n 
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with the flux factor 

F = (2ir) 5 yj fa ■ p*)* - m*. (32) 

The factors 1/4 and emerge for the simple reason that the polarization of initial and 

final particles is not considered. Since the relative phases among various meson exchanges in 
the amplitude of Eq. (26) are unknown, the interference terms are ignored in our concrete 
calculations. 

In the amplitude given above, we have neglected the pp initial-state-interaction (ISI) 
factor which has been used by many authors in investigating pp — > ppr] reaction [42, 43]. The 
role of the ISI in the production of a heavy meson from an NN collision is basically to reduce the 
cross section near threshold by an overall factor with little energy dependence, although with 
some controversy about its absolute value. While the ISI factor with pp — > ppr\ reaction near 
threshold is relatively simple with 3 P pp initial state dominant, it becomes technically more 
complicated for pp — > NNn with several important initial partial waves even at energies very 
close to the threshold. In practice, for the pion production reactions, the ISI dumping effect is 
not included explicitly in previous calculations [22, 23] and is done by adjusting the vertex form 
factors to reproduce the total cross sections. Since how to treat this dumping factor does not 
influence our main conclusion here, we follow the simple practical approach as these previous 
calculations [22, 23] on the same reaction, although the ISI could be taken into account by 
some more complicated approaches, such as by including the partial wave ISI factors [42, 43] 
or distortions of wave functions [44, 45, 46]. 
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3 Numerical results and discussion 

With the formalism and ingredients discussed in the former sections, we evaluated the total 
cross section versus the kinetic energy of the proton beam (T P ) for the pp — > pnn + reaction by 
using the code FOWL from the CERN program library, which is a program for Monte Carlo 
multi-particle phase space integration weighted by the amplitude squared. The results for Tp 
ranging from 0.2 to 1.3 GeV are shown in Fig. 2 and Fig. 3 along with experimental data [47, 48] 
for comparison. 

In Fig. 2, the results obtained from with and without including the 3 S\ np FSI factor 
are represented by solid and dashed curves, respectively. We can see our theoretical result 
without FSI agrees well with the experimental data at beam energies near 1.0 GeV. However, 
at lower beam energies,the calculated total cross sections underestimate the data by a factor 
of 4 or more. In view of the important role of FSI for the near-threshold enhancement, we 
also tried to include it and found that the result with FSI is indeed in excellent agreement 
with the experimental data over a wide range of beam energies. Both the energy dependence 
and the absolute magnitudes of the experimental data are very well reproduced. Note that 
in reality, the np 3 S\ partial wave should dominate only at energies close to threshold. For 
energies of a few hundred MeV above threshold, the contributions from higher partial waves 
become important. The FSI effect due to the strong np 3 Si interaction should become less 
important as energy increases. Hence after taking this effect into account, one would expect 
to get theoretical results closer to the solid curve at low energies and to the dashed curve at 
higher energies in the Fig. 2. 

In Fig. 3, individual contributions corresponding to A(1232) with tt + and 7r° exchange, 
A r *(1440) with a and 7r° exchange, and nucleon pole with a and 7r° exchange are shown in 
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comparison with the experimental data by dotted, short-dotted, dot-dashed, dot-dot-dashed, 
dashed, and short-dashed curves, respectively. The contribution from the A ++ (1232) pro- 
duction by the ir + exchange is found to be dominant for T P between 0.6 and 1.1 GeV. The 
A + (1232) with it exchange gives one order of magnitude smaller contribution for the cross 
section than A ++ (1232) with n + exchange over the whole energy region, simply due to smaller 
isospin factor of 1/9. In Ref. [23], by only including the contribution of the A(1232) resonance 
the authors reproduced the experimental points for present channel fairly well, however, visibly 
underestimated the data points close to threshold. Actually, from Fig. 3, one can see that 
the nucleon pole by a meson exchange plays a more important role for beam energies close to 
threshold. This is consistent with the previous calculation for the pp — > ppir reaction [45], 
where it is shown that the prediction with only pion exchange [46] underestimates the data by 
a factor about 5. And there are more examples from other reactions showing the significance 
of the sub-threshold nucleon pole contribution, such as in J/ip — > pnir + [17, 49]. In contrast, 
contribution of the nucleon pole by n exchange is much smaller than by a exchange over the 
whole energy region. Furthermore, it is worth noting that contribution of the iV*(1440) reso- 
nance by a exchange becomes quite significant at beam energies above 1.1 GeV, larger than all 
other contributions except A ++ (1232) contribution, though it is negligibly small at lower beam 
energies. The contribution of the iV*(1440) by ir° exchange is much smaller by a factor of about 
10 or more, which indicates that by no means can it account for the clear Roper peak in the 
invariant mass M n7T + spectrum. Anyway, the a-meson exchange is found to play a significant 
role, which has also been pointed out in a previous study of ap reaction [18]. Possible contri- 
bution from the neighborhood iV*(1535) resonance is also checked and is found to be negligible 
in the present energy region. 

In our present calculation, we have not included the contribution from nN S-wave re- 
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scattering process, which was found to be the dominant production mechanism for the 7r + 
production at the energies very close to the threshold [50]. Since the main interests of the 
present work is for the N* production at higher energies and in our approach the uncertainties 
from FSI and ISI are also quite large at the energies very close to the threshold, we refer the 
readers who are interested in the very near threshold regime to the more dedicated study of 
this regime by Ref. [50]. 

In our concrete calculation, the interference terms between various resonances are ignored 
because the relative phases from various resonances are unknown. In order to get some feeling 
about the importance of the interference terms, we calculated the magnitudes of a few largest 
ones. At the high energy end T P = 1.3 GeV, the largest two contributions come from A ++ 
and iV*(1440) resonances as shown in Fig. 3. For the cross section, the contribution of the 
interference term from these two resonances is found to be in the range of ±5% of the contri- 
bution from the A ++ resonance alone. At the lower energies around 0.4 GeV, the largest two 
contributions from A ++ and off-shell nucleon are comparable, but the contribution from their 
interference term is found to be only around ±1% to the total cross section by integrating over 
the whole three-body phase space. 

The invariant mass spectra and Dalitz plot for the final particles of pp — > pmr + reaction 
at T P = 1.3 GeV are also calculated as shown in Fig.4. While the A ++ (1232) peak dominates 
the pir + invariant mass spectrum, a peak due to iV*(1440) is clearly visible in the n7i + invariant 
mass spectrum. As a rough comparison with data, preliminary results with two different triggers 
and without proper acceptance correction from CELSIUS-WASA Collaboration [20, 21] are also 
shown in the pn + and nn + invariant mass spectra. Individual contributions to the invariant 
mass spectra from A(1232), iV*(1440) and nucleon pole are also given in Fig. 4 by dashed, 
dotted, and dot-dashed curves, respectively. From the nn + invariant mass spectrum, one can 
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see that the reflection of A ++ (1232) resonance gives a large but flat contribution, the isospin 
suppressed A + (1232) contribution gives a weak peak around 1.23 GeV, and the Roper iV*(1440) 
resonance is absolutely needed to reproduce the peak around 1.36 GeV. In the calculation for 
T P = 1.3 GeV, the n-p FSI factor is not included since the data for both prr + and nn + mass 
spectra have no evidence for the FSI effect which should cause further enhancement for both 
mass spectra at high energy end. The phenomena may suggest that at beam energy about 1.3 
GeV the contribution from p-n higher partial waves has already exceeded the 3 Si partial wave 
to be the dominant contribution and the S-wave FSI becomes unimportant. The forthcoming 
acceptance corrected results from the CELSIUS- WASA Collaboration and future experiments 
at COSY with the newly installed WASA-at-COSY detector [51] or at HIRFL-CSR with the 
scheduled in hadron detector [52] will be very helpful for constraining the ingredients in our 
model calculations. 

In our calculation, the only free parameter we adjust is the for the AnN vertex 
form factor. Previous calculations used = 0.65 GeV [23] and 0.63 GeV [53]. Whereas they 
included only the A resonance to reproduce the total cross sections for the pp — > pnn + reaction, 
we need to reduce it to 0.59 GeV to allow contributions from the off-shell nucleon pole and the 
Roper resonance for the best description of the total cross section. 

The experimental study [20, 21] and our theoretical study here suggest that the pp — > 
pnn + reaction provides a very good place for studying the N* resonances. Since all the contri- 
butions from A* + resonances are well constrained by the corresponding A* ++ contributions by 
an isospin scaling factor of 1/9 and hence suppressed, the mr + invariant mass spectrum for this 
reaction provides a rather good isospin filter for N* resonant peaks. Therefore the study of this 
reaction should be extended to higher energies at COSY and HIRFL-CSR with 4ir detectors to 
provide various differential cross sections and the Dalitz plot. It will definitely help us to look 
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for those "missing" A* ++ resonances with large coupling to p + p and iV* resonances with large 
coupling to Ncr or Nuj at higher energies. The information on various couplings of observed A* 
and N* resonances may also help us gain some insight on the nature of these resonances [32]. 
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Appendix 



In this appendix we present explicitly the expression of each amplitude in our calculation of 
Eq.(26) as the following: 



M(A(1232),tt h 



fit- 



>l(A(1232),7r ) 



M(iV*(1440),7r ) 



M(N*(14A0),cr) 



M(p(93S),n°) 



M{p(93B),a) 



G n (k n )u n {p n , s n )7 5 j£ n u 2 (p2, s 2 ) - (exchange term with p x <-> p 2 ), (33) 

= ^^9l N ^ N {kl)F^{kl)u n {p n , s n )plG^\q)KuM, *i) 
3 p 

G n (k n )u 3 {p 3 , s 3 )75 hu 2 {p2, s 2 ) - (exchange term with p x <-> p 2 ), (34) 

= ^^^ JV /r(^)^ ,JV (^)«n( P n, S)! )75 ^( 144 °)(g) 75 JUlCPl,*) 
771 tt 

G^(K)u 3 (p 3 , s 3 )7 5 jt w u 2 (p2, s 2 ) - (exchange term with p x <-> p 2 ), (35) 
= v^^^^7v^.^(^)Fj v ^(^)« n ( Pn , Sn ) 75 ^G^( 144 °)(g)«ibi, si) 
G a (k a )u 3 (p 3 , s 3 )u 2 (p 2 , s 2 ) - (exchange term with p x <-> p 2 ), (36) 

G 7T (k 7T )u 3 (p 3 , s 3 )7 5 ^U2(p2, s 2 ) - (exchange term with pi <-> p 2 ), (37) 
= v / 2^^(F^(^)) 2 F(g)« n (p n , Sn )75 ^G N {q) Ul {p u 8l ) 
G a (k a )u 3 (p 3 , s 3 )u 2 (p 2 , s 2 ) - (exchange term with pi <-> p 2 ) (38) 



with u n (p n , s n ), u 3 (p 3 , S3), «i(pi, Si), w 2 (p 2 , s 2 ) the spin wave functions of the outgoing neutron, 
proton in the final state and two initial protons, respectively, p^ and Av(o-) stand for the 4- 
momenta of the outgoing n- meson and the exchanged 7r(<r)-meson, respectively, pi and p 2 
represent the 4-momenta of the two initial protons. 

The amplitudes corresponding to the "pre-emission" graphs are not given here, but are 
included in actual calculations. With these amplitudes, the J2 Si S S/ l-^| 2 m Eq.(31) is computed 
with the standard MATHEMATICA package. 
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Table 1: Relevant paramters of A(1232) and JV*(1440). 

Figure 1: Feynman diagrams for pp — > pnir + reaction. Note that there are also the pre- 
emission counter parts of these diagrams where pion is emitted before collision. These graphs 
are not shown here, but included in the calculations. 

Figure 2: Total cross section vs T P for the pp — > pmr + reaction from our calculation 
compared with data [47, 48]. The solid and dashed lines represent the results with and without 
3 S\ np FSI, respectively. 

Figure 3: Contributions of various components and their simple incoherent sum as a 
function of T P for the pp — > pnn + reaction compared with data [47, 48] , without including np 
FSI. The dotted, short-dotted, dot-dashed, dot-dot-dashed, dashed, and short-dashed curves 
represent contributions from A(1232)(7T + and 7r° exchange), iV*(1440)(er and tt° exchange) 
and nucleon pole (a and 7r° exchange), respectively; the solid curve represents their simple 
incoherent sum. 

Figure 4: Invariant mass spectra and Dalitz plot for the pp — > pnn + reaction at T P = 
1.3 GeV, compared with the preliminary data with two different triggers and without proper 
acceptance correction (open circles from Ref. [20] and solid circles from Ref.[21] with arbitrary 
normalization). The dashed, dotted and dot-dashed lines stand for individual contributions 
from A(1232), iV*(1440) and nucleon pole, respectively, while the solid line represents their 
simple incoherent sum. In the Dalitz plot, the size of the squares indicates the magnitude of 
the number of events. 



23 



Table 1: 



Resonances Width(MeV) Decay Channel Branching ratios(%) g 2 

A*(1232) 118 vriV 100 19.54 

iV*(1440) 300 ttN 65 0.51 

aN 7.5 2.84 
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Figure 1: 
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Figure 2: 
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■ □□□□□□□□□□□□ = 
■ □□□□□□□□□□□□□□□ ■ 
= □□□□□□□□□□□□□□□□□ ■ 
□ □□□□□□□□□□□□□□□□□□□ ■ 
. □□□□□□□□□□□□□□□□□□□□□ ■ 
■ □□□□□□□□□□□□□□□□□□□□□□□ 
□□□□□□□□□□□□□□□□□□□□□□□□ = 
= □□□□□□□□□□□□□□□□□□□□□□□ □ □ 
■ □□□□□□□□□□□□□□□□□□□□□□□□□□ ■ 
□ □□□□□□□□□□□□□□□□□□□□□□□□ □□ □ 
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o □□□□□□□□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□□□□□□□□□□□□□□□□□□□□□□□□□□□ - 
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